Claims : 

1. A method of processing a digital image , comprising the 
steps of: 

providing digital data indexed to represent positions on 
a display, said digital data being indicative of an intensity 
value I i (x,y) for each position (x,y) in each i-th spectral 
band; 

adjusting said intensity value for said each position in 
each i-th spectral band to generate an adjusted intensity 
value for said each position in each i-th spectral band in 
accordance with 

N 

Y, w n ( lo g IiiX'V) - log[I ± (x,y) *F n (x,y)]) , i = l,...,S 

23 = 1 

where S is the number of unique spectral bands included in 
said digital data and, for each n, W n is a weighting factor 
and F (x,y) is a unique surround function applied to said each 
position (x,y) and N is the total number of unique surround 
functions ; 

filtering said adjusted intensity value for said each 
position of said image in each of said S spectral bands using 
a filter function wherein a filtered intensity value R n -(x,y) 
is defined; and 

selecting a maximum intensity value V^x^y) from the 
group consisting of said intensity value I^XjY) and said 
filtered intensity value R-(x,y). 
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2. A method according to claim 1 wherein each said unique 
surround function is a Gaussian function. 



1 3. A method according to claim 2 wherein said Gaussian 

2 function is of the form 

-r 2 

e ^ 

3 satisfying the relationship 

-r 2 

k n f f e c?n dxdy = 1 

jjSt where 

|j r = sjx 2 +y 2 

'*5 and, for each n, k n is a normalization constant and c n is a 

ff 6 unique constant for each of said N unique surround functions. 
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4. A method according to claim 1 further comprising the step 
of multiplying said filtered intensity value R,-(x,y) by 



log 



Bl^x.y) 
E Ji(x # y) 



to define a color-restored intensity value R',(x,y), where B 
is a constant, wherein said step of selecting using said 
filtered intensity value R,(x,y) is replaced with the step of 
selecting a maximum intensity value V^x^) from the group 
consisting of said intensity value l,-(x,y) and said color- 
restored intensity value R' r (x,y). 



5. A method according to claim 1 wherein said each position 
(x,y) defines a pixel of said display. 

6. A method according to claim 1 wherein, for each n, W n =l/N. 

7. A method according to claim 1 further comprising the step 
of displaying an improved image using said maximum intensity 
value V, (x,y) . 

8. A method according to claim 4 further comprising the step 
of displaying an improved image using said maximum intensity 
value V,- (x,y) . 



9. A method according to claim 1 further comprising the steps 
of: 

defining a classification of said image based on dynamic 
range of said image in each of said S spectral bands; and 

selecting said filter function based on said 
classification of said image. 

10. A method according to claim 9 wherein said step of 
defining comprises the step of using image statistics 
associated with said image in each of said S spectral bands. 

11. A method according to claim 10 wherein said image 
statistics include brightness and contrast of said image in 
each of said S spectral bands. 

12. A method according to claim 4 further comprising the 
steps of: 

defining a classification of said image based on dynamic 
range of said image in each of said S spectral bands; and 

selecting said filter function based on said 
classification of said image. 

13. A method according to claim 12 wherein said step of 
defining comprises the step of using image statistics 
associated with said image in each of said S spectral bands. 



1 14. A method according to claim 13 wherein said image 

2 statistics include brightness and contrast of said image in 

3 each of said S spectral bands. 
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15. A method of processing a digital image, comprising the 
steps of: 

providing digital data indexed to represent the positions 
of a plurality of pixels of a J-row by K-column display, said 
digital data being indicative of an intensity value I(x,y) for 
each of said plurality of pixels where x is an index of a 
position in the J-th row of 'said display and y is an index of 
a position in the K-th column of said display wherein a JxK 

image is defined; 

convolving said digital data associated with each of said 
plurality of pixels with a function 

c 2 

e c 

to form a discrete convolution value for each of said 
plurality of pixels, said function satisfying the relationship 

kjfe ° 2 dxdy = 1 

where 

r = yjx 2 +y 2 

k is a normalization constant and c is a constant; 

converting, for each of said plurality of pixels, said 
discrete convolution value into the logarithm domain; 

converting, for each of said plurality of pixels, said 
intensity value into the logarithm domain; 



20 subtracting, for each of said plurality of pixels, said 

21 discrete convolution value so-converted into the logarithm 

22 domain from said intensity value so-converted into the 

23 logarithm domain, wherein an adjusted intensity value is 

24 generated for each of said plurality of pixels; 

25 filtering said adjusted intensity value for each of said 

26 plurality of pixels with a filter function wherein a filtered 

27 intensity value R(x,y) is defined; and 

28 selecting, for each of said plurality of pixels, a 
|| maximum intensity value V(x,y) from the group consisting of 
H said intensity value I(x,y) and said filtered intensity value 
H R(x,y) . 

0\ 16. A method according to claim 15 wherein the value of said 

C| constant c is selected to be in the range of approximately 

fie 0.01 to approximately 0.5 of the larger of J and K. 

' i 17. A method according to claim 15 further comprising the 

2 step of displaying an improved image using said maximum 

3 intensity value V(x,y). 

1 18. A method according to claim 15 wherein said step of 

2 filtering includes the step of selecting said filter function 

3 based on dynamic range of said JxK image. 
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1 19. A method of processing a digital image, comprising the 

2 steps of: 

3 providing digital data indexed to represent the positions 

4 of a plurality of pixels of an J-row by K-column display, said 

5 digital data being indicative of an intensity value I-(x,y) 

6 for each i-th spectral band of S spectral bands for each of 

7 said plurality of pixels where x is an index of a position in 

8 the J-th row of said display and y is an index of a position 

9 in the K-th column of said display wherein a ( JxK) i image is 
10 defined for each of said S spectral bands and a JxK image is 
H defined across all of said S spectral bands; 

Jj convolving said digital data associated with each of said 

|j| plurality of pixels in each i-th spectral band with a function 

Co 



5| for n=2 to N to form N convolution values for each of said 

|S plurality of pixels in each said i-th spectral band, said 

if 

m function satisfying the relationship 

-i 2 

k n J J e c?n dxdy = 1 



17 where 

r = yjx 2 +y 2 

18 and, for each n, k n is a normalization constant and c n is a 

19 unique constant; 

20 converting, for each of said plurality of pixels in each 
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said i-th spectral band, each of said N convolution values 



22 into the logarithm domain; 

23 converting, for each of said plurality of pixels in each 

24 said i-th spectral band, said intensity value into the 
2 5 logarithm domain; 



subtracting, for each of said plurality of pixels in each 
said i-th spectral band, each of said N convolution values so- 
converted into the logarithm domain from said intensity value 
so-converted into the logarithm domain, wherein an adjusted 
intensity value is generated for each of said plurality of 
pixels in each said i-th spectral band based on each of said 



|2 N convolution values; 

f| forming a weighted sum for each of said plurality of 

f| pixels in each said i-th spectral band using said adjusted 

§fe intensity values; 

|£ filtering said weighted sum for each of said plurality of 

i|7 pixels in each said i-th spectral band with a filter function 

1*1:1: 

38 wherein a filtered intensity value Ri(x,y) is defined; and 

39 selecting a maximum intensity value V,-(x,y) from the 

40 group consisting of said intensity value I,(x,y) and said 

41 filtered intensity value R,-(x,y). 

1 20. A method according to claim 19 wherein the value for each 

2 said unique constant c n is selected to be in the range of 

3 approximately 0.01 to approximately 0 . 5 of the larger of J and 

4 K. 
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21. A method according to claim 19 further comprising the 
step of multiplying said filtered intensity value R f (x,y) by 



[ Bl^x.y) 
log 

3 to define a color-restored intensity value R',-(x,y), where B 

4 is a constant and S is a whole number greater than or egual to 

5 2, wherein said step of selecting using said filtered 

6 intensity value R,(x,y) is replaced with the step of selecting 

7 a maximum intensity value V,- (x,y) from the group consisting of 
r| said intensity value I,-(x,y) and said color-restored intensity 
0 value R ' (x, y) . 

f| 22. A method according to • claim 19 further comprising the 

8^ step of displaying an improved image using said maximum 

is 

% intensity value V^XjY). 
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1 23. A method according to claim 21 further comprising the 

2 step of displaying an improved image using said maximum 

3 intensity value V-(x,y). 

1 24. A method according to claim 19 further comprising the 

2 steps of: 

3 defining a classification of said JxK image based on 

4 dynamic range of each said (JxK) ^ image; and 

5 selecting said filter function based on said 
^ classification of said JxK image. 

% 25. A method according to claim 24 wherein said step of 

B| defining comprises the step of using image statistics 

% associated with each said (JxK) - image. 

m 

Hi 2 6. A method according to claim 2 5 wherein said image 

g| statistics include brightness and contrast of each said (JxK) i 

3 image . 

1 27. A method according to claim 21 further comprising the 

2 steps of: 

3 defining a classification of said JxK image based on 

4 dynamic range of each said ( JxK) 1 - image; and 

5 selecting said filter function based on said 

6 classification of said JxK image. 
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1 28. A method according to claim 27 wherein said step of 

2 defining comprises the step of using image statistics 

3 associated with each said ( JxK) f image. 

1 29. A method according to claim 28 wherein said image 

2 statistics include brightness and contrast of each said (JxK) , 

3 image . 



30. A method of processing a digital image, comprising the 
steps of: 

providing digital data indexed to represent positions of 
an image having S spectral bands for simultaneous output on a 
display, said digital data being indicative of an intensity 
value I,(x,y) for each position (x,y) in each i-th spectral 
band; 

adjusting said intensity value for said each position in 
each i-th spectral band to generate an adjusted intensity 
value for said each position in each i-th spectral band in 
accordance with 

N 

Y,W n dog Ii(x,y) - log [I^x.y) *F n (x,y) ]) , i=l,...,S 

n=l 

where S is a whole number greater than or equal to 2 and 
defines the total number of spectral bands included in said 
digital data and, for each n, W n is a weighting factor and 
F n ( x f y) i s a unique surround function of the form 



e 



satisfying the relationship 

-z 2 

k n f je c " dxdy = 1 



where 



r = \jx 2 +y 2 

and, for each n, k n is a normalization constant and c n is a 



19 unique constant where N is the total number of unique surround 

20 functions; 

21 filtering said adjusted intensity value for said each 

22 position in each i-th spectral band with a filter function 

23 wherein a filtered intensity value R,-(x,y) is defined; 

24 multiplying said filtered intensity value R,(x,y) by 

log [ B /' U ' y) ' 

_ 2 = 1 

25 to define a color-restored intensity value R^-fxj), where B 
2;6 is a constant; and 

|| selecting a maximum intensity value V-(x,y) from the 

|| group consisting of said intensity value I^x.y) and said 

f§ color-restored intensity value R'^x^) • 

(1% 

Q 3i. a method according to claim 3 0 wherein, for each n, 

f| W n =l/N. 

" \ 32. a method according to claim 30 wherein the value for each 

2 said unique constant c n is selected to be in the range of 

3 approximately 0.01 to approximately 0.5 of the larger of J and 

4 K. 

1 33. A method according to claim 30 further comprising the 

2 step of displaying an improved image using said maximum 

3 intensity value V ; (x,y). 
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1 34. A method according to claim 30 further comprising the 

2 steps of: 

3 defining a classification of said image based on dynamic 

4 range of said image in each of said S spectral bands; and 

5 selecting said filter function based on said 

6 classification of said image. 

1 35. A method according to claim 34 wherein said step of 

2 defining comprises the step of using image statistics 
Q associated with said image in each of said S spectral bands. 

Qt 36. A method according to claim 3 5 wherein said image 

L- :; 2 statistics include brightness and contrast of said image in 

*3 each of said S spectral bands. 
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